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Abstract 

We establish a correspondence between the resummation of world lines and the 
diagonalization of the Hamiltonian for a strongly correlated electronic system. For 
this purpose, we analyze the functional integrals for the partition function and 
the correlation functions invoking a slave boson representation in the radial gauge. 
We show in the spinless case that the Green's function of the physical electron 
and the projected Green's function of the pseudofermion coincide. Correlation and 
Green's functions in the spinful case involve a complex entanglement of the world 
lines which, however, can be obtained through a strikingly simple extension of the 
spinless scheme. As a toy model we investigate the two-site cluster of the single 
impurity Anderson model which yields analytical results. All expectation values 
and dynamical correlation functions are obtained from the exact calculation of the 
relevant functional integrals. The hole density, the hole auto-correlation function 
and the Green's function are computed, and a comparison between spinless and 
spin 1/2 systems provides insight into the role of the radial slave boson field. In 
particular, the exact expectation value of the radial slave boson field is finite in 
both cases, and it is not related to a Bose condensate. 

Key words: Strongly correlated electrons, Field theories, Non-perturbative 
methods 

PACS: 71.27.+a, ll.10.-z, 11.15.Tk 



* Corresponding author 

Email address: Raymond.Fresard@ensicaen.fr (Raymond Fresard). 
1 Present address: LASMEA, UMR CNRS-Universite Blaise Pascal 6602, 24 avenue 
des Landais, 63177 Aubiere Cedex, France. 



Preprint submitted to Nuclear Physics B 



6 February 2008 



1 Introduction 



In contemporary solid state research, strongly correlated electrons comprise 
the most fascinating albeit intangible physical systems. They cover a wide 
range of phenomena, including high temperature superconductivity, colossal 
magnetoresistance, aspects of the fractional quantum Hall effect, and even 
electronic reconstruction in oxide electronic devices which are built on in- 
terfaces of strongly correlated films. Whereas their importance is generally 
perceived, a fundamental comprehension is still not achieved, especially for 
high-temperature superconductivity. 

This unfortunate absence of a well established theoretical scheme or even solu- 
tion is not surprising: strong electronic correlations are based on (sufficiently 
strong) local interactions in real space but the Fermi surface, the concept on 
which the physics of metals is firmly rooted, is defined and understood in mo- 
mentum space. Correspondingly, a theoretical investigation is either built on 
momentum or real space approaches which allow to treat either the kinetic 
(band) term or the interaction accurately. However, a momentum-space weak 
coupling approach is insufficient to generate the desired new energy scales and 
fix points whereas standard perturbation theory from the highly degenerate 
local (atomic) limit suffers from severe drawbacks [l|. 



Nevertheless, many effective strong coupling theories expand, in a generalized 
sense, with respect to local models. A model with a single local interaction 
term is the Anderson impurity model. It is the prominent strong-coupling 
many-body model which can still be solved exactly (with certain restrictions) 
and which has been understood in basically all aspects (see, for example 
0; 0; 0; 0; 0|)- It may justly be seen as the paradigm of a strongly correlated 
many-body system. A successful scheme to investigate lattice models with 
on-site interactions originates from a self-consistent extension of the Ander- 
son model. The self-consistency is generated through a dynamical mean-field 
theory (DMFT) which singles out a site with strong local interaction; this 
site couples to an electronic bath, the effective medium, the local density of 
states of which is calculated self-consistently @;S]. Actually, the DMFT is ex- 
act for infinite space dimensions, a limit which was introduced in Ref. [9| for 
correlated electron systems. However, it is missing the spatial correlations. In 
recent years it has been devised to treat clusters which can couple to various 
bath systems in order to investigate correlations with a spatial extension of 



the cluster size 10; 11; 12; 13; 14 



In our theoretical study we will focus on a different approach, the slave boson 
technique 15|; [l6|; [17| . The formalism entails a local decomposition of electronic 
excitations into charge and spin components. Electron creation and annihila- 
tion operators are thereby represented by composite operators which separate 



2 



into canonical operators with bosonic and fermionic character. However these 
operators are enslaved in the sense that their respective number operators 
have to fulfill a local constraint. The original idea was to decouple spin and 
charge degrees of freedom; other, modified schemes attribute to each type of 
excitation a bosonic mode which allows to study the correlated system in a 
saddle point approximation 



17; 18; 19 



This mean field approach has been 
successful when set against numerical simulations: ground state energies [2(| 
and charge structure factors show excellent agreement 2l|, as the procedure 



is exact in the large degeneracy limit [18j; |22 



While the saddle point approximation allows to calculate translationally in- 
variant expectation values in momentum space, the corresponding mean field 
solution is not a priori legitimate. The objection is concerned with the local 
decomposition of the electron field into fermion and slave boson components. 
This implies that the model acquires a local gauge invariance with the conse- 
quence that Elitzur's theorem [23J prevents the (slave) bosonic fields to acquire 
a non-zero expectation value. In fact, it is the phase fluctuations of the boson 
field which suppress the finite value or condensation of these fields. 



One alternative to avoid such a condensation has been devised by Kroha, 



Wolfle, and coworkers [2J; |25(. In their approach the local gauge invariance 
is guaranteed through Ward identities in a conserving approximation. The 
projection onto the physical sector of the Fock space is achieved with an 
Abrikosov procedure by sending the Lagrange multiplier of the constraint 
to infinity] 2 1 The other alternative is to use a radial decomposition of the 



bosonic field [27j, the details of which were presented in a previous paper by 



two of the authors 28]. In the limit of large on-site interaction, the bosonic 
fields in radial representation reduce to their respective (real) amplitude as the 
time derivatives of the conjugate phase can be absorbed in a time-dependent 
Lagrange multiplier field. 



In this article we provide a scheme for the solution of cluster models in radial 
slave-boson representation. We present in sufficient detail the calculation of 
correlation and Green's functions for a two-site cluster of the single impurity 
Anderson model, in order to exemplify our scheme. Although the model can 
be diagonalized without slave boson technique we esteem the explicit solution 
in the radial decomposition of considerable significance. First, it relates the 
world line expansion of slave boson path integrals to the quantum states in 
the Fock space, in particular for entangled states. This is achieved through a 
decomposition of the fermionic determinant into resolvents at each time step. 
Second, it allows to compare these exact results (e. g., for the slave boson 



amplitude) to saddle point evaluations and to assess their validity 29 



Note that a similar procedure can be set up without introducing slave bosons [261 ] . 
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The article is organized as follows: in Sec. [2], we introduce the functional in- 
tegral formulation of the two-site cluster model. We give expressions for the 
action, partition function and hole density as well as for the hole autocorre- 
lation function in terms of radial slave bosons. The spinless system is studied 
first in Sec. [3] where, through the derivation of the partition function, we show 
how to proceed from a world line representation to the representation with 
quantum states in the Fock space. In Sec. HI we show how our formalism allows 
to derive results for the spinful case from a straightforward extension of the 
spinless case. The Green's function necessitates a distinct derivation of the 
fundamental connection between the slave-boson path integral representation 
and the hamiltonian scheme which is the object of Sec. [51 



2 Functional integral formulation of the two-site cluster model 

2.1 Hamiltonian and radial slave boson representation 



The single impurity Anderson model (SIAM) has been investigated with a 
variety of techniques and for many different purposes. One of them consists of 
testing a new approach, in particular against exact results. Here we adopt a 
similar spirit in order to link the evaluation of the path integral representation 
of thermodynamic and dynamical quantities to their computation through a 
straightforward diagonalization of the Hamiltonian. For the SIAM it reads: 



^ = Ei^+e c )c fci(7 +E4^+^E( c iA + h - c -) +f/ II 44* 

k,cr a k,a cr=1\i 



where U is the on-site repulsion, which is hereafter taken as infinite. The 
operators c\ a (c k a ) and d\ (d a ) describe the creation (annihilation) of the 
band electrons and impurity electrons respectively, with spin cr; the kinetic 
energy in the band is denoted t k , and e c and are the band and impurity 
energy levels, respectively. The hybridization is given by V. 

The link between the two schemes is actually a complex procedure when the 
impurity is coupled to an infinite bath. In order to lower this complexity 
we reduce the bath to a single site, in which case the diagonalization of the 
Hamiltonian is easy and all relevant results can be obtained analytically. Nev- 
ertheless, the problem is non-trivial when handled in the functional integral 
formalism. The level of difficulty depends on the functional representation 
which is used. For our purpose, a promising one is that of the slave boson 



(SB) representation in the radial gauge [28(. It is based on the original repre 
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sentation by Barnes [l5( and is augmented in that respect that the underlying 
U(l) gauge symmetry, originally discussed by Read and Newns 271 ]. is fully 
implemented, as the phase of the bosonic field is integrated out from the out- 
set. Accordingly, the original field d a is represented in terms of a real and a 
Grassmann field for each spin: 



d n ,cr X n+lfn,cr 
yi,cr ay nJn,cn 



(2) 
(3) 



where x n and x n+1 are the slave boson field amplitudes at time steps n and 
n + 1, and f n a is the auxiliary fermion field. The shift of one time step for x in 
the relation for d n a is necessary to obtain a non-zero value of the Grassmann 
integration for the Green's functions — (d (V)dt(O)) as clearly shown for its 
calculation in the atomic limit in Ref. [28). More precisely, the path integral 
is zero if x n+l is replaced by x n in Eq. fl2]). Moreover, Eqs. fl2]) and (j3J) are 
required in order to properly represent the hybridization term in the action 
as given below. Further detail on this matter can be found in Ref. 28 . 



2.2 Action and partition function 



Following Ref. [28J, the path integral representation of the partition function 
of the two-site cluster is given by: 



z = j im fn / UDiujijDic^cij r r d Xn ) e ~ s (4) 



where the action S may be written as the sum of a fermionic part, Sf, and a 
bosonic part, with 



S f = J2 S L° = J2 [ C n, CT (Cn,a - ^cC„-l, CT ) + f^M^a ~ ^nfn-\,a) (5) 
a n=l cr 

+ V5x n {c ] n) J n -i^ + fl^Cn^) 

N 

Sb = J2 I 5 ( iX n(x n - 1) + Wx n (x n - 1))] , 



n=l 



where L c = e 5< - ec >*\ L n = e s( - €d ^+* A ") = L d e lSXn , A n is the time-dependent 
real constraint field, n denotes the time steps, and 5 = /3/N, with (3 = l/k B T 
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and N the number of time steps. Here, Sf (Sf,a) is bilinear in the fermionic 
fields, and the corresponding matrix of the coefficients will be denoted as 
[S] ([So]). The positive real number W is sent to infinity at the end of the 
calculation which guarantees the projection onto the physical subspace. The 
above treatment of the bosonic field is specific to radial slave bosons: no phase 
variable appears, and the above form cannot be obtained by transformations 
of the conventional functional integral in the Cartesian gauge without further 



assumptions (see Ref. |28| ) 



Inserting Eq. (jSj) into Eq. (TjJ yields a particularly suggestive expression of the 
partition function Z in the functional integral formulation, of the kind: 



Z= lim V\...Vn det [S] , (6) 



where detfS] is the determinant of the fermionic matrix defined by Eq. (jHJ) 
V n is defined as: 



Vn= I^S^L r°°dx n e -*[<A»(*n-i)+w*n(x„-i)] j (7 ) 

oo 27T J— oo 



and acts as a projector from the enlarged Fock space "spanned" by the aux- 
iliary fermionic fields down to the physical one. Explicitly, the action of these 
projectors on the various contributions resulting from det [S] are found to be: 



V n -l = l, (8) 

"Pn ' J'n 1) (9) 

V n -L n = L d , (10) 

V n -L n x n = 0, (11) 

V n -L 2 n = 0, (12) 

V n -x 2 n = l. (13) 



As will be seen below no further property of V n will be needed for our pur- 
pose. We note that there is some freedom in writing the projectors V n , and 



alternative expressions exist 28(. However the properties Eqs. (15)1— (IT31) are 



independent of the particular form of V n . 
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3 Application to the spinless fermion case 



We first consider a spinless fermion system for simplicity. Even though this is 
a non-interacting problem, the level of complexity of its path integral repre- 
sentation following from Eqs. (jl])-©, is equivalent to the one of a fully inter- 
acting problem. The matrix representation of the action S/ jCr of such system 
is a 2N x 2N square matrix whose explicit expression in the basis {c n>a , / ni<T } 
reads: 



[So] 



I 2 

[£ 2 ] I 2 



(14) 



[£ N ] l 2 / 



where I2 is the 2x2 identity matrix and [C n ] are 2x2 blocks given by: 



-L c 5Vx n ^ 
5Vx n —L n ) 



(15) 



at time step n, 1 < n < N. Note that the matrix [So] as defined in Eq. ([T 
has the same structure as the action matrix S^ in Chap. 2 of Ref. [30]. 



3. 1 Partition function 



The partition function Zq of the spinless fermion system has a form similar 
to that of Eq. (jSJ) except that det [S] is replaced by det [So]. Its calculation is 
straightforward since we only have to evaluate: 



Zq = hm V\ 

N-+00 

W^oo 



2N 

V N £sgn(P)nSp(n),n 
{P} n=l 



(16) 



where sgn(P) is the signum function of permutations P in the permutation 
group S2N and Sjj are the matrix elements of [So]. Since the time step n is 
only involved in Sj j2n -3 and S ii2 n-2 we m ay recast Eq. (|T5|) into: 
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Z = lim Vsgn(P) V 2 (Sp(i),iSp(2),2 

W->oo {P} 



X 



V N 



At this point it is straightforward to verify that performing the projections 
only implies to make use of Eqs. (|8l)-( |T0i) . We are left with: 



2N 



Z ° = i im Yl s S n ( p ) II S P(n),n = lim det [S' ] , 



{P} 



n=l 



N^oo 



where SI • are the elements of the 2N x 2N matrix [S' ] defined as: 



[Sol 



2 



-[£] 



(19) 



In Eq. (fT9~|) the 2x2 matrix blocks [C] are similar to the blocks [C n ] except 
that L n becomes L d , and x n is replaced by 1: 



[4 



5V -L H 



(20) 



In the form of Eq. (1181) it is now obvious that Zq can be readily obtained. 
We notice that Eq. fll9p is the expected action matrix for this free fermionic 
problem. Besides, it is straightforward to extend the above calculation to the 
case of an arbitrary bath. Unfortunately the calculation becomes considerably 
more involved in the spin 1/2 case, which leads us to develop another strategy 
to that purpose. We first present it in the spinless case, before extending it to 
the spinful case. 



3.1.1 Generation and resummation of the world lines 

Part of the difficulty in computing Zq is that the time steps are mixed in 
the fermionic determinant, in contrast to the bosonic part of the action repre- 
sented by the projectors V n , Eq. ([H]). Therefore, transforming this determinant 
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into a form where the time steps are decoupled, is desirable. Achieving this 
amounts to handle all the world lines following from the action in Eq. (151) 
which represents a problem equivalent to a particle in a time- dependent field 
with a time-dependent hopping amplitude. 

In order to generate the dynamics of the world lines we first expand det [So] 
along the first two columns. We obtain: 



det [So] = 1 x M 2 2 + L C M 3 2 2 + 5Vx 2 Ml 2 - SVx 2 Ml 3 - L 2 M\ A (21) 
+ (L C L 2 - (5Vx 2 ) 2 ) M 2 4 . 



Here the notation is as follows: we construct a matrix similar to [So], but we 
only include time steps 1 and m > n. M"- is a minor of this matrix, where 
both the z-th and j-th rows, together with the first and second columns, are 
removed. 

At this stage we may proceed with the generation of the world lines. To that 
aim we need to express the minors M 2 as linear combinations of the minors 
M 3 which in turn can also be expressed as similar linear combinations of the 
minors M 4 , and so forth up to the time step N. The recurrence relation that 
we have established takes the following form: 



(m? :2 \ 

Ml 2 

Mf A 
\ M lJ 



"#(n+l),l" 



'^(n+l), 



where the four matrix blocks 



'K^ n) ' 1 ] =(1), 

j^(n),2l _ |"^(n),3 

~K^]=(L c L n ), 



L c 5Vx r , 
5Vx n L n 



(23) 
(24) 
(25) 
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at time step n, 1 < n < N, define the 6x6 block diagonal matrix K^] . GO 



The matrices 



describe the evolution of the two-site system along the 
world lines at each time step n. Iterating this procedure up to time step N 
yields the determinant det[So] in the following scalar product form: 



det [S ] 



1, L r , 5Vx 2 , -5Vx 2 , -L 2 , L C L 2 



N 



(«) 



n=3 



1 

5V Xl 
-SVx x 
-U 
\ L C L X J 



(26) 



where the row vector is identified from Eq. (|2~IT) and the column vector has 
been obtained from the last time step of the iteration process. Since the minus 
signs in Eq. fliZEl) cancel, they can be discarded for further considerations. 



The above expression corresponds to the full resummation of the world lines, 
some of them being represented in Fig. [TJ The first contribution, labeled by 
1, corresponds to the subspace with zero electron while the last one, labeled 
by L c L d , corresponds to the subspace with two electrons. In both cases the 
world lines are "straight": namely no hopping process takes place, and the 
system is right away in an eigenstate. In Eq. (1261) they correspond to the 



terms involving 



and 



respectively. The structure of the world 



lines in the one-electron subspace is more intricate. 

In order to gain an intuitive picture of these world lines let us first consider the 
trivial functional integral representation of an interactionless electron where 
the /-field directly represents the physical electron. We begin with the pro- 
cesses where the electron is on the "band" site at time step one. If it stays there 
during all time steps, the resulting contributions Z^ to Z will be given by 
= Un=i L c -> namely there is one factor L c per time step and the world line 
is straight. If, on the contrary, the electron hops onto the impurity at time step 
m, and back to the "band" at time step m', the corresponding contribution 

z£l, c to Z results in Z^, c = (lC=i ^o) SV (ut^+i L A ) 5V (utm'+i L ) • 
Higher order processes in V follow accordingly. Complementary processes are 
those where the electron resides on the impurity at time step one. The con- 
tributions of all these processes to Z will be weighted by both the number of 



In Eq. (|25p a term of order 8 2 , that vanishes in the limit N — > oo, is neglected. 



10 



ooooooooooo 




Fig. 1. Typical world lines of a two-site system. Thick (thin) lines denote the el- 
ementary processes with the electron sitting on the impurity ("band") site. Full 
(empty) circles denote the elementary processes with a hole sitting on the "band" 
(impurity) site. Dashed lines represent hopping processes. The horizontal axis runs 
with the time steps, while in the vertical direction different sites are displayed. 

hopping processes and the difference in energy between the two levels. Assum- 
ing €d < e c results in L c < and for world lines involving the same number 
of hopping process, the world line containing the largest number of factors 
yields the largest contribution. 



If we now return to our representation the contributions of the world lines to 
the partition function follow in a similar fashion, except for that i) the factor 
corresponding to hopping process at time step m is given by SVx m , and ii) 
the particle line at time step n corresponding to the electron sitting on the 
impurity site results in a factor L n . Accordingly the contribution of the world 
line labeled by L c in Fig. [T] to det [So] is: 

/m-l \ / m'-l \ / N \ 

n L c \5vx m n m wx m > n m > (27) 

\n=l / \n=m+l J \n=m'+l / 

while the world line labeled by in Fig. [1] yields 



/m— 1 \ / m' — l \ 

Y[L n \5Vx m i I] L c \5Vx r 

\n=l / \n=m+l / 




(2f 



3.1.2 Connecting with the hamiltonian language 

We observe that the structure of Eq. (|26|) is not manifestly translationally 
invariant in time. For this reason, we proceed to bring Eq. (|26l) into a form 
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where no particular time step is singled out. We note that the (row) column 
vector in Eq. (1261) can be identified with the (rows) columns of the matrix 
blocks - disregarding the cancelling minus signs. Accordingly, we can 

rewrite Eq. (1261) as: 



det [S ] = 


E 

{a2,...,a N } 




0203 


• ' a N _ 1 a N ""-ajvl 


+ 


E 

{a2 r --,Q!iv} 




/^(3),2 

«2«3 ' 


^(JV),2 ^(1),2 
• "ajv-idff ajvl 


+ 


E 

{a 2 ,...,a N } 


/r (2),3 
n 2a 2 


Zf(3),3 
± ^a 2 a 3 ' 


^(1)>3 


+ 


E 

{a 2 ,...,a N } 


^la 2 


K (3),4 
- tl a2a3 


■ " [i ajv„ia]v "ivl 



since the matrices 



Kin) 



are block diagonal and symmetric. Observe that the 



first and last lower index in line three of Eq. (129]) is not 1 but 2. 



In Eq. (1291) . the first sum is equal to 1; the second and third sums are the 

/ 



diagonal elements of the matrix product rinLi 



L r 5Vx r 



5Vx n L n 



respectively. 



The last sum is equal to Yln=i L n - Therefore Eq. (|29|) reduces to the trace 
of a 4x4 matrix that is the product of the N block diagonal matrices [JC n ] 



whose elements are 



and 



respectively: 



N 



det [So] = Tr J] [JC„\ , 



n=l 



where the matrix [K. n ] is given by: 



(30) 



[lC n ] 



L r 5VX T 



SVx n L n 



L c L n 



(31) 



Here, it is of interest to note that we have established a direct link between 
the world line picture embodied in the 6x6 matrices K^ 1 ' and the simpler 
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description in terms of quantum states in the Fock space through the 4x4 ma- 
trices [JC n ]. Indeed, when performing the world line expansion using Eq. ( j26|) . 



the propagation matrices 



acquire an involved structure following from 
the initial conditions attached to the world lines. This is best seen in the one- 
electron sub-space: at time-step one, the electron may either be on site c or on 
site d (see Fig. [1]), and the corresponding dynamics is governed by the matrices 



and 



contribution to det 



. Summing up all these world lines yields the one-electron 
Sq]. Thus, one needs to handle two 2x2 matrices, one in 
each case. In contrast, a single 2x2 matrix needs to be treated on the level of 
Eq. (J3DJ). 



In the above form of the fermionic determinant, Eq. (|30|) . the time steps are 
decoupled, which greatly simplifies the projection onto the physical Hilbert 
space. Indeed, in terms of the matrix [q] given by: 



[q]=V n ([JC n ]) 



L c 5V 
5V L d 



L C L, 



(32) 



we obtain the partition function Zq as: 



N 



Zn = lim Tr TT \q] = lim Tr 

TV^oo J " L 1 11 N^oo 
n=l 



1a- 8 







e c - fx -V 
-V e d - /I 



e c + e d - 2/i j 



N 



(33) 



Namely, we recover here the hamiltonian matrix in the Fock space. There- 
fore the expansion of det [Sq] in minors together with the recurrence relations, 
Eq. (12"2"|) . allows for a correspondence between the ensemble of the world lines 
and the hamiltonian matrix. It is apparent that the complexity of this inter- 
relation depends on whether or not the system is in an eigenstate at time 
step one. If this is the case, there is one single straight world line, and the 
connection is obvious. Otherwise there is a proliferation of world lines, here 



controlled by the matrices 



and 



which recombine to yield the 



contribution resulting from the 2x2 block of the matrix [fC n ]. A result equiv- 
alent to Eq. (133|) was already obtained by Barnes [3lj], though in a totally 
different fashion. 
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3.2 Hole density 



In our path integral formalism, the expectation value of the amplitude of the 
slave boson field at time step m is (x m ). It simply represents the hole density 
1 — rid which can be written as (x m ) = (6^6 m _i) where b represents the original 
Barnes slave boson, and (x m ) is finite. In contrast the expectation values of 
the boson operators are zero: (6^) = (b m ) = 0, for each time step because of 
the fluctuations of their respective phase factor, in agreement with Elitzur's 
theorem. Note that expectation values of higher order moments of also 
non zero: (x^) = ( x m) 7^ for any real positive parameter a. 

In the case of spinless fermions we calculate (x m ) as: 



Z Q {x m ) = lim Vi . . . V N (det [So] x m ) (34) 



lim^ V 1 ...V N I x m Ti J[ [K n ] J 



If we introduce the 4x4 matrix [Qx] = V n {x n [JC n ]) for all n, and [Q] the 4x4 
diagonal matrix that satisfies [q] = [Uq] [Q] [Uq]\ [Uq] being the eigenvector 
matrix, Eq. (134"]) becomes: 



Z (x m ) = lim Tr ([Q]^ 1 [U Q f [Q x ] [U Q ]) . (35) 



Using Eqs. (I8l)- ffl~0l) . and the 5^0 limit, we obtain that the matrix [Qx] 
reduces to the representation of the hole density operator on the impurity in 
the Fock space: 



[Qx]ij = ki SjA + 5 ^ (36) 

Correspondingly, Eq. fl33|) expresses (x m ) as the averaged value of the hole 
density operator on the impurity, represented by its matrix elements in the 
basis of the eigenstates of the Hamiltonian ([C/g] [Qx] [Uq]) and weighted by 
the Boltzmann factors [Q]^ 1 . Therefore, in contrast to a Bose condensate, 
(x m ) is generically finite and may only vanish for zero hole concentration. Its 



numerical evaluation will be presented in a forthcoming paper 29 
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3.3 Density- density correlation function 



To obtain further insight into the approach, it is of interest to compute dy- 
namical correlation functions. Probably, the simplest one is provided by the 
hole density autocorrelation function on the impurity site. When expressed in 
terms of eigenstates of the Hamiltonian, it takes the form: 



((1-mimS)) (l-n d (5))) = -1 ^ e s( - m ~ N)Ea e~ mSE -''\(ip a \l - n d |Vv)| 2 , 

u a,a' 

(37) 

where the eigenvalues E a can be obtained from the diagonalization of the 
Hamiltonian. The latter can be read from, e. g., Eq. ([3"3"j) . 

In our path integral representation, the autocorrelation function (xix m ) may 
be written as: 



Z (xix m )= lim Pi... Pat (det [So] xix m ) (38) 

N— >oo 



lim Vx ... V N XiX m Tr TT [/C„] 

AT— >oo \ 1 / 



which reduces to 



Zo(x 1 x m ) = lim Tr ([Q x ] [q] m - 2 [Q x ] [qf-™) , (39) 

after application of the projectors, Eqs. fl8|)- f|T3|) . If we now introduce the 
eigenvalues and eigenvectors of the Hamiltonian, Eq. fl39|) can be recast into: 

Z ( Xl x m ) = Jim Tr ([Q] N ~ m [U Q f [Q x ] [U Q ] [Q] m - 2 [U Q ? [Q x ] [U Q ]) . (40) 



In this form we recognize the standard expression in Eq. (l37|) . Indeed, the 
combination [Uq]' [Q x ] [Uq] represents the matrix elements of the hole density 
operator in the basis of the eigenstates of the Hamiltonian, and the factors 
[Q] the exponential factors in Eq. f l37l) . 
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4 Spin 1/2 system 



We turn now to the spinful case. In this section, we show how results obtained 
for the spinless system are relevant and useful to derive in a straightforward 
fashion the corresponding quantities of the spin 1/2 system. 



4-1 Partition function 



Since up- and down-spins are decoupled, and since the fermionic contribu- 
tion to the action is identical for both of them, we can write the fermionic 
determinant as: detfS] = det[So] 2 . This allows to express the determinant as 
the product of the traces of the matrix products II^Li [K-n] • By making use of 
the mixed product property of the Kronecker product: {[A] [C]) <g> ([£>] [D]) = 
([A] ® [B]) ([C] ® [D]), we obtain: 



N 

det [S] = Tr J] [JC n ] ® [JC n ] . (41) 

71=1 

The partition function is obtained by combining Eqs. (1411) and Since the 
time steps are decoupled, the evaluation of the integrals over the x and A fields 
is straightforward. The tensorial products [JC n ] ® [JC n ] yield 16x16 matrices, 
which, after application of the projectors V\ . . . Vn, become block diagonal 
with a 4 x 4 zero block. They all reduce to the same 12x12 real symmetric 
matrix [k]. The matrix [k] represents the projection of the tensor product 
[JC n ] <S> [JC n ] for all n as 



[% = V n {[K, n \ iuh [/C n ] i2i . 2 ) (42) 



with the convention 



'4(n - 1) + 1 2 if (1 < %\ < 2 and 1 < i 2 < 4) 

and (i! = 3 and 1 < i 2 < 2), (43) 
= i 2 + 10 if (n = 4 and 1 < i 2 < 2). 

and similarly for j. The remaining matrix elements of the 16 x 16 matrix form 
a vanishing 4x4 separate block and have been discarded in Eq. (1421) . Finally 
we obtain a simple expression for the partition function Z of the two-site 
single impurity Anderson model: 



Z = lim Tr [k] N . (44) 

N—>oo 



16 



with 



L c 5V 
5V L d 



\ 



L c Ld 



L c 




^ 
Ll L r 5V L r 5V 



L C 5V L c L d 




8V 
L C 5V 





L 2 C L A 











L d 
L c L d 





















L c L d 
L\L d ) 



(45) 



When expanded to lowest order in 5, the blocks of [k] represent the hamiltonian 
matrix in the Fock space, in the same fashion as in Eq. (|33|) . Diagonalizing 
these blocks yields the expected expression of the partition function Z: 



Z = 1 + 3 exp [-0 (e c + e d - 2/x)] + 2 exp [-0 (2e c + e d - 3/x)] 



+ 2 ex P 

+ exp 
3 -=±i 



e c + ed 



- H + j< 



e d \ 



+ V 2 



+ 2^ 2 



(46) 



Through the exact calculation of the functional integrals we have recovered 
this result that can also be derived from the diagonalization of the Hamilto- 
nian. Note that V is multiplied by different coefficients in the one-particle and 
two-particle states. In contrast to the spinless case Eq. ( 1331) the eigenvalues 
of the hamiltonian matrix entering Eq. ( )45l) result from entangled states, the 
entanglement being achieved by the projection onto the physical Fock space in 
Eq. (142]) . Note that we here obtained the hamiltonian matrix without having 
explicitely used any basis of the Fock space. It naturally arose as the projected 
tensor product of the basis appropriate to the spinless case. 
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4-2 Hole density and autocorrelation function 



In the spinful case, the hole density (x m ) is given by: 



Z(x m ) = lim V x . . . V N (det [S\ x m ) (47) 

N— »oo 

= lirn^ V 1 ...V N ( x m Ti J[ [K n ] <%> [JC n ] J 



For the spin 1/2 case, the counterpart of the matrix [Qx] is the matrix 
[ICx] = V n (x n [K. n ] ® [/C n ]). It is a 12x12 matrix, the elements of which may 
be expressed as: 



[fcxlij = + $i,2 $j,2 + $i,5 83,5 + 5j,B, (48) 

using Eqs. (!Hl)- (|T3l) . in the 5 — > limit. In this form it represents the hole 
density operator on the impurity site in the Fock space. Thus Eq. (T4T1) becomes 



Z(x m )= hm Tr {[K x ] [k] N ~ l ) , (49) 

which has exactly the same form as for the spinless case and a similar inter- 
pretation applies. 

As for autocorrelation functions such as (xix m ), we can adopt the same pro- 
cedure to obtain 



Z(x 1 x m ) = lim Tr ([Kx] M"^ l^x] [k] N ~ m ) , (50) 

which, when again introducing the eigenstates of the Hamiltonian, can also 
be identified to the ordinary expression in Eq. (1571) . On the formal level, the 
use of the slave boson representation in the radial gauge greatly simplifies the 
evaluation of the dynamical correlation functions of the operators that can be 
represented by the radial slave bosons. 



5 Green's function 



We turn to the impurity Green's function G a (q—p). When expressed in terms 
of the eigenvalues of the Hamiltonian, E a — which can be read from, e. g., 
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Eq. (H| for spin 1/2 



- it reads: 



ZGM-V) =-£ e^- Ar -*^<Va|da|V'«'>e* (p - ff)B -'<^|4^«> > (51) 



In the radial gauge the creation and annihilation operators are expressed in 
terms of auxiliary fields as given in Eqs. (j2J) and (j3J). We obtain G a (q — p) as 
follows: 



N 



ZGM - p) = - lim J] / II ^[/n^. /Ul^K^, (52) 

/"°° 6d\n f°° , _ 5 , ff 

O aX n e Jq,(rJp,cj X q+l X pi 

J — oo ^7T J — oo 



in the language of functional integrals. Note that the three other Green's func- 
tions involving the three expectation values: (d s ), {c d) a ) and (c cj ) , 
can be calculated in the same fashion, except that they are simpler to evaluate 
since they contain at most one amplitude of the slave bosonic field, x, unlike 
the one we chose to study. 



5.1 Derivation of the pseudofermion Green's function 



Following standard procedures (see, e. g., Negele and Orland (30|), we cast 
Eq. ( 152]) into the form: 



ZG a {q -p) = - lim Vx...V N (det [So] G P , g x p x q+1 ) , (53) 

N — >oo 

where Q VA is the minor of one of the matrix elements of the 2x2 block that 
shares the same row labels as [C p ] and column labels as in the matrix 

[Sq] as defined in Eq. (fUjl (if q — N then the block to be considered is [£i]). 
The minor Q VA is the unprojected pseudofermion Green's function. 

For the subsequent calculations we set p = N — m + 1 and q = N. The minor 
Gn-tu+in can be calculated as: 



d 

Q N - m +i,N = q- det [S(a; m)] , (54) 



where [S(a; m)] is given by: 
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[S(a;m)] i:j = [So] id + aS ii2 (N-m+i)Sj,2N- 



(55) 



To calculate det [S(a; m)], we find it convenient to move the last two columns 
to the left in the matrix [S(a; m)\ in Eq. (1551) . In the same fashion as for det [So], 
we expand det [S(a; m)] along the first two columns. Once the derivative of 
det [S(a; m)] with respect to a is calculated we obtain the three contributions: 



d_ 

da 



det [S(a; m)] = -1 x M} a + L c M\ a + 5V Xl ML, 



(56) 



where the three minors M.\ a , M\ c(L and M\ a are defined in the appendix. 
They can be expressed as linear combinations of M.\ a , Ml a and Mj a in 
the same fashion as for the minors defined in the previous sections. They also 
follow a recurrence relation that reads: 



Ml c>a 

\ M la ) 



1 
L, 
SVx 



n+l 





5Vx n+1 
L, 



\ ( M n+1 \ 



J n+1 J 



mW 



(57) 



up to time step n = N — m. Here we recognize the blocks defined in Eqs. 
and f!24p entering the matrix [i^ n+1 )] which we encountered during the eval- 

^tf and Mf^ are linear 



uation of the partition function. Then, M.i a 
combinations of the minors 



N ~ m+1 and Mz A m+l : 



f M N-m\ 



M N-m 



l 8,a 



-1 
L, 
5Vx 



N-m+l 



M N-m+l 



M N-m+l 



(5f 



l 3,4 



The minors M.\^ and M.^^ are tightly related to M" 4 and M^ A : they are also 
built on a matrix similar to [So], but with the difference that only the time 
steps m with m > n are included. M^i ^ s a mrn o r of this matrix, where the first 
two columns and both the i-th and j-th rows are removed. Accordingly, the 
recurrence relations for Aii^ m+1 and ■M^ m+l are also given in Eq. f[2"2"j) and 
hence we also need to introduce the minor Aii^ m+1 . Thus, their evaluation 
involves the blocks of Eqs. (|24|) and ([25]) . which enter the matrix [i^ < - n+1 ' ) ] 
jointly with the last three components of the column vector in Eq. (|26|) . taken 
at time step N. Therefore, to calculate GN-m+i,N as defined above we have 
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to consider the following set of six minors: Aii ja , -M-L c ,a, -Ms,a for time steps 
1 < n < N — m, and and for time steps N — m + 2<n<N. 

Combining the above steps, we can write the unprojected pseudofermion 
Green's function as: 



'JV-m+lJV 



Tr 



K 



(TV) 



X 



K 



(N-m+l) 
P 




(59) 



where the time steps enter in decreasing order. The 6x6 matrices 



K 



(N-m+l) 
P 



K 



(N) 



representing the annihilation (creation) of a fermion in the 



world line language are given by 



K 



(N) 
f 



K 



(N-m+l) 
P 



1,3 



'■J 



SVx N 5 iy i 5j A + L N 5 iy2 5j£ + L C L N S it2 5 jfi 
■ 1 Si,5 + L c S ifi 5 j>2 + 5Vx N - m+ i 5 ifi S jj3 . 



(60) 



In Eq. (J59J), 



is a 6 x 6 block diagonal matrix the non zero elements of 



which are the two blocks 



The matrix K 



and 



i^(n),2l defined in Eqs. (|23l) and (j22j). 



is also 6x6 block diagonal matrix determined by 



K 



(n) 



K 



(n) 



(61) 



Equation flBTJl) can also be interpreted on the basis of Fig.[TJ In fact, computing 
(d N a d N _ m+l CT ) can be visualized as the resummation of particular subsets of 
world lines. They are naturally split into sets involving for iV — m + 2 < 

is controlling the dynamics of the world lines; 



n < N — 1 , in which case 



K 



(n) 



and sets excluding for 1 < n < N — m, in which case 



K 



{n) 



controls the 



dynamics. The transition from the first (second) subset to the second (first) 
is taken care of by the matrix m+1 ' ) ( ) . 



This expression for the Green's function in the world line language can also 
be related to its counterpart in the hamiltonian language. Indeed, the trace of 
the matrix product in Eq. (1591) can be written in terms of 4x4 matrices: 
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Q N -m + l,N = Tr [F N \ X [] l^n ) >< [^N-m+l] X \ ]J [K* 

\ \n=N-l J \n=N-m ) ) 

(62) 

where 



m+lj 



^oo o o^ 


10 

\0 L c 5Vx N _ m+1 J 



and [Fn] = 



( 5Vx N L N \ 
L C L N 




(63) 



/ 



characterize the creation and annihilation of an electron, respectively. The 
matrices [/C<] and [/C>] result from [/C n ] as: 



(64) 



5.2 Spinless case 



Now that the unprojected pseudofermion Green's function has been converted 
to a compact form, we can evaluate the physical Green's function which, in 
the spinless case, reads: 



Z G (q - p) = - lim Vi...V N Tr 

W-*oo 



N-m+2 
XiX N - m+ i [F N ] X ( Yl 

n=JV-l 



X [<t>JV- m+ l] X J] M 



\n=N —m 



(65) 



With the application of the projectors Pi . . . Pjv, we obtain: 



ZoG (q -p) = - hm Tr ([Q x ] [F] [q>] m ' 2 [0] [q^' 1 ) , (66) 
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where [Qx] is given by Eq. (|36|) . The matrices [0] and [J 7 ] are defined as 
[0] = V N - m+ i (x N _ m+1 [$7v_ m+ i]) and [J 7 ] = V N {[F N ]) respectively. They can 
be read off from Eq. (I6"3l if x is replaced by 1, and L n by L d . Note that 
[(f)] = [J 7 ] only in the limit 5^0 when L c — > 1 and L<j — > 1. In the limit 
5^0, they coincide with the matrix representations of the operators f* and 
/ respectively. The matrices [q K ] and [g > ] are given by [g < ] = V n ([JC<]) and 
= respectively, and are easily related to [q]: 



=[«]«- ^i^i- (67) 

we can now reshape Eq. fl66l) into the form: 

ZoG (q -P) = - lim Tr \[Q<] N - m ^ ([U Q ]* [Q x ] [J 7 ] [U Q ]) (68) 

x [Q>] m - 2 f[C/ Q ]t [0] [Uq] 



where [Q < ] and [Q > ] are obtained from the diagonalization of [q K ] and [g > ], re- 
spectively. Note that the same matrix [Uq] diagonalizes the matrices [q] , [q K ] , 
and [g > ]. 

In Eq. ( |68|) . the product [Uq]^ [Qx] [F] [Uq] is the representation of the anni- 
hilation operator d in the basis of the eigenstates of the Hamiltonian, while 
[Uq]^ [0] [Uq] represents d\ as can be easily verified explicitly in the limit 
5 — > 0. The factors [Q > ] and [Q < ] determine the time evolution. Therefore 
Eq. fl68|) can be easily identified with Eq. fl5T|) . Note that all eigenvalues seem 
to contribute to the Green's function in Eq. fl5Tl) . and only the matrix ele- 
ments of the creation (annihilation) operators restrict the set of eigenvalues 
that effectively contribute to the Green's function. In contrast, in Eq. (|68|) 
some of these restrictions are contained in the factors [Q > ] and [Q < ] which 
replace the full set of eigenvalues, that would be contained in the matrix [Q]. 

It is tempting to compare the physical electron Green's function (including 
the factors x) to the projected pseudo-fermion Green's function (without the 
factors x). Straightforward algebra yields: 



lim Vn-m+l (xN-m+l [$jV-m+l]) = lim V N - m+1 ([$iV-m+l]) , (69) 
o— >0 o^O 



and 
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lim V\Pn 

<5^0 



(*i [K<] [F N ]) 



lim V{P N {[K<] [F N ]) 



(70) 



Therefore, as a particularity of the spinless case, the factors x in Eq. 0651) 
play no role, and both Green's functions coincide. Consequently the same 
result for the physical electron Green's function would have been obtained 
by substituting x by 1 in Eqs. (jSJ) and (j3J), and accordingly in the fermionic 
contribution to the action Sf. Incidentally, such a procedure is in complete 
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agreement with the original suggestion by Kotliar and Ruckenstein 
modify the expression of the physical electron operator by introducing square 
root factors, when extended to the spinless case. 



5.3 Spin 1/2 system 



Again, as shown below, results obtained for the spinless system can be im- 
mediately applied to derive the Green's function in the spinful case. Inserting 
Eq. d59]) into Eq. ([53]) yields: 



ZG a (q - p) = - lim Vi...V N XiX N ^ m+ i 

iv— >oo 
W-*oo 



(N-m+2 
II [£„]®[/C; 
n=N-l 

X ([K N -m+l] ® [$N-m+l]) X ( II i K n] ® [K 

\n=N —m 



(71) 



which after application of the projectors V\. . . Vn becomes: 



ZG a {q -p) = - lim Tr ([K x ] [£] [k>} m ~ 2 [<p] [k<] N - m ~ l ) , (72) 

N—>oo > ' 

with the matrices [ip] = V N - m+1 (x N _ m+1 [/C;v-m+i] ® [§ N - m +i}) and [f] = 
"Pn ([JCn] ® [-Pjv])- Leaving out entries which do not contribute in the limit 
iV — > oo they read: 



Mij = 1 <5i,3 S jt i + L c (5 iA 5 jt2 + S i:7 5 j:5 ) + L\ 5 ii8 6 ji6 

[Qij = L c 8i tl S j>3 + L c L d (5 ij2 b~ jA + 5 i;5 Sjj) + L 2 c L d S ifi 5 jt8 . (73) 

The matrices [fc < ] and are given by [A; < ] = V n ([)C n } Cg> [/C<]) and [A; > ] = 

T^ndJCn] ® [/C^]) respectively. The asymmetry in the representation of the 
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physical electron creation and annihilation operators (Eqs. and (j3J)) is also 
apparent in Eq. (172!) . Indeed the operator d) a is represented by the matrix [ip] , 
and d a by the product of the matrices [JCx] [£]■ In contrast to the spinless case 
the factors x in Eq. (I7T|) play a role, and the projected pseudofermion Green's 
function and physical electron Green's function differ. 



6 Conclusion 



In summary we have established a new scheme which provides a fundamental 
connection between the representation of expectation values and dynamical 
correlation functions in the hamiltonian language and their counterpart in the 
slave-boson path integral formulation. This has been achieved for the U = oo 
spin- 1/2 single impurity Anderson model through their exact evaluation for a 
two site cluster. The new scheme allowed us to compute the partition function 
and the hole density, expressed as the expectation value of the radial slave 
field x. Moreover the Green's function and the hole auto-correlation function 
were evaluated within this scheme. 

We verified that the exact expectation value of the slave boson amplitude 
field x is finite, as postulated in mean-field calculations, even in this extreme 
quantum case. It is therefore not related to the condensation of a boson, 
which would necessarily vanish in such a calculation. The suppression of the 
condensation originates in the use of the radial representation, where the phase 
of the boson is integrated out in the first place. We note that higher slave boson 
correlation functions such as (x n (r)x m (0)) reduce to (x(t)x(O)). Therefore 
the field x bears little resemblance to ordinary complex bosonic fields. The 
corresponding calculations follow a similar scheme as those for the partition 
function. 



Through an independent calculation we obtained both the physical electron 
and pseudofermion Green's functions. In the spinless case, the projected pseud- 
ofermion Green's function is finite, and it is identical to the one of the physical 
electron. Therefore a "perturbation theory" -like factorization of the latter as 
a product of the boson and pseudofermion Green's functions does not ap- 
pear appropriate in general, but it may still be valid in particular frequency 
ranges, such as the low frequency domain. In the latter case, a mean-field 
decoupling looks more appropriate. It is likely to provide a better agreement 
with the exact result if the square root factors, originally proposed in [17| , are 
introduced. 



It is also of great importance to understand that our formalism allows imme- 
diate and straightforward use of results, which were obtained for the spinless 
system, in order to derive those of the spinful case: the proposed scheme first 
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treats the coherent states of fermions in the two spin sectors (up and down 
spin) separately. The world lines of particles with different spin projection 
evolve independently. Only in a final step, when the full fermionic determi- 
nant is built from the product of the determinants of the two spin species, 
the projection onto the physical space with no double occupancies is straight- 
forwardly achieved through the projection rules, Eqs. fl8l- ffT3l . applied to the 
entries of the determinant. Here, the projection is easily accomplished in the 
Fock space which, when directly enforced for the world lines, produces a com- 
plicated entanglement of coherent states. For larger systems the exact resum- 
mation of the world lines and their respective projection — as presented in, 
e.g., Eqs. (I39I) . (!50|) and (166]) — is difficult on the analytical level, but probably 
not on the numerical level. As an alternative to the exact calculation one may 
consider a plain saddle point approximation scheme when tackling spatially 
extended systems. Unfortunately the latter fails to reproduce the exact result 
even for the two-site problem, and it will be necessary to determine appropri- 
ate quasiparticle weight factors from the two-site solution within an effective 



slave boson approach, similar to the Kotliar-Ruckenstein scheme [17J. This 



challenge will be addressed in future work [34 



The extension of the above scheme to the spin rotation invariant formulation 
of the t-J model where the phases of all the bosonic fields can be gauged away 
is desirable [18]. Work along this line is in progress. 

One may wish to extend such a calculation to other representations of this 



model, such as the ones based on Hubbard X-operators [32|. This unfortu- 
nately poses another challenge, since the "angular part" of the respective 
action is intrinsically off diagonal in time, which makes the integral over 
the angular variables significantly more difficult. This also holds true for the 
Kotliar-Ruckenstein representation where one of the bosonic fields is complex 



18l ; 1171 ; [33j. Alternatively one may also consider weak-coupling approaches, 
such as the Hubbard-Stratanovich decoupling of the interaction term in the 
charge channel. Even if it were possible to evaluate the partition function ex- 
actly, using the corresponding form of Eqs. ([6]), (J7J) and fT4"TT) . it would still 
require a major effort to obtain dynamical response functions. Nevertheless 
such a calculation deserves further study. 
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A Expression of the minors M.\ a , •M.'i a and M.g a 

The minors M. \ a , M. \ c a and A4j a in Eq. (|56|) are obtained from the expansion 
of det [S(a; m)\ given below for p = N — m + 1. Their definition follows as: 

• M\ a by expanding det [S(a; m)] along the two first columns and eliminating 
the '(2N - l)th and 2pth lines. 

• M.\ a by expanding det [S(a; m)] along the two first columns and eliminating 
the second and 2pth lines. 

• M\ c a by expanding det [S(a; m)} along the two first columns and eliminat- 
ing the first and 2pth lines. 

In the expansion of det [S(a;m)} M.\ a is multiplied by [S( y a;m)] 2N _ 11 = 

1 and — [S(a; m)] 2 2 = —a, M.\ a is multiplied by —[S(a;m)] 21 = SVxi 

and [S(a;m)} 2 2 = a, and M-\^ a is multiplied by [S(a; m)] 1 1 = L c and 
[S(a;m)} 2 2 = a. The three minors above satisfy a recurrence relation given 
in Eqs. ( IBTl) and ( J58l) . while det [S(a]m)} reads: 



det[S(a;m)]= (A.l) 
L c -5Vxt 1 
-5Vx 1 Lt 1 
— L c 5V%2 
5Vx 2 —L 2 










a 



-L c SVx p 



5VX T 



1 






1 



-L c 8Vx p+ i 
5Vx p+ i —L p+ i 



1 





1 



1 
1 

— L c 5Vxn 
SVxn —Ln 
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